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Abstract 


A  cross  flow  drag  coefficient  model  for  finite  length  body-strake  eombinations  has  been 
developed.  The  basis  for  the  new  model  is  transonie  test  results  of  two-dimensional  body- 
strake  eombination  by  Maeha.  Adjustment  for  finite  body  was  done  based  on  the 
empirieal  eorrelation  for  eireular  bodies  by  Jorgensen.  Additional  adjustment  eonsiders 
the  relatively  large  eross  flow  drag  eoeffieient  of  shake  alone  at  subsonie  Maeh  numbers. 

The  model  was  ineorporated  in  eomponent  buildup  predietion  method.  The  stability 
derivatives  were  estimated  using  the  Missile  Dateom  eode.  The  inviseid  eontributions  to 
the  normal-foree  and  pitehing-moment  eoeffieients  are  based  on  these  stability 
derivatives.  The  contribution  of  the  separated  eross  flow  eonsiders  body  sections  with 
shakes  as  one  unit  and  uses  the  new  drag  eoeffieient.  The  method  was  applied  to  seven 
eonfigurations  at  Maeh  numbers  between  0.9  and  6.8.  In  most  oases  good  agreements  was 
obtained  between  prediotions  by  the  Missile  Dateom  eode,  the  new  method  and  test  data. 
The  present  method  gives  better  results  than  the  Missile  Dateom  code  for  two 
eonfigurations  that  feature  reotangular,  or  mostly  reotangular,  shakes.  It  was  observes 
that  when  the  eross  flow  Maeh  number  is  higher  than  1.0,  the  present  method 
overestimates  the  normal-foree  and  the  stability  margin.  This  implies  that  the  present 
eross  flow  drag  eoeffieient  model  needs  further  adjustment. 
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Nomenclature 


a  defined  by  Eqn.  (B.6.a) 

A  defined  by  Eqn.  (B.6.b) 

b  maximum  inelusive  span 

B  width  of  eonfiguration 

e  ehord  of  strake 

Cdc  eross-flow  drag  eoeffieient  of  a  finite  body 

Cdo  eross-flow  drag  eoeffieient  of  cireular  eylinder  (2-D) 

Cdo  eross-flow  drag  eoeffieient  of  finite  cireular  body 
Cd|  cross-flow  drag  coefficient  of  isolated  strake  (2-D) 

Cdi  cross-flow  drag  coefficients  for  s/d^l  .0  (2-D) 

Cd2  cross-flow  drag  coefficients  for  s/d=2.0  (2-D) 

Cd(|,  cross-flow  drag  coefficient  of  body-strake  cross  section  (2-D) 

Cd(|)  cross-flow  drag  coefficient  of  finite  length  body-strake  combination. 

Cl  lift  coefficient 

Cm  pitching-moment  coefficient 

Cn  normal-force  coefficient 

d  body  diameter 

dR  reference  diameter,  maximum  body  diameter 

f  length  of  body 

M  Mach  number 

s  inclusive  span  of  strake 

Sp  plan-form  area 

Sr  reference  area,  (7i:/4)dR 

X  longitudinal  coordinate 

Xo  location  of  tip  of  strake 

xr  reference  point  for  pitching-moment  and  center  of  pressure 
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Greek 


a  angle  of  attack 
Subscripts 

1  contribution  of  inviscid  flow 

2  contribution  of  separated  cross-flow 

New.  Newtonian  theory 

Designation  of  the  components 

B  body  alone 

B-S  body-strake  combination 

S  strake  alone 

SU  strake  unit,  strakes  with  mutual  effects  with  the  body 
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Introduction 


In  recent  years  there  is  a  growing  interest  in  missile  configurations  that  feature  strakes, 
rather  than  wings.  The  advantage  of  strakes  is  their  small  spans,  resulting  small  volume, 
which  is  beneficial  for  internal  carriage  and  launching  from  canisters. 

Codes  that  are  based  on  component  buildup  methodologies  estimate  the  linear  and 
nonlinear  characteristics  of  the  components  separately.  Thus,  the  nonlinear  contributions 
of  bodies  and  wings  are  estimated  using  different  methods  or  databases.  Flow 
visualization  of  body-strake  combinations,  by  Verle,^  shows  that  the  leeward  vortex  wake 
of  body-strake  combination  is  common  to  the  two  components,  as  depicted  in  Fig.  1.  This 
implies  that  the  contribution  of  the  separated  cross  flow  to  the  normal-force  of  the 
combination  may  be  different  than  the  sum  of  the  contributions  of  the  components. 

Literature  survey  yielded  only  one  experimental  study  of  the  cross-flow  over  straked 
circular  cylinders.  Macha^’  ^  measured  the  cross-flow  drag  coefficient  of  plain  cylinders 
and  in  combination  with  strakes  at  Mach  numbers  between  0.6  and  1.0  and  span  to 
diameter  ratios  up  to  2.0. 

The  objective  of  the  present  study  is  to  explore  the  benefit  of  using  unified  body-strake 
cross-flow  drag  coefficient  in  the  estimation  of  the  normal-force  and  center  of  pressure 
location  of  such  configurations. 


t»cfW8  IwoStej*  #  ■Ufou/  d0  tu5«/<o«s  cylinanquts  comportant  das  nagaoiraa  mlncas  parmuas  ou  das  moua- 

-  Cyllndfa  i  net  goihiqua  poiniu  sana  (abc)  ai  avac  nagaoiraa  <a'bc') . 

®  I'axtradoa  (aa')  at  tranchaa  iransvaraafaa  aiiuOas 

=  UOO  i  mi-nagaoira  (bb'  r  =  12  D)  at  t  rafri4ra  (cc'  x'  »  18  D) 


Fig.  1:  Flow  visualization  of  the  leeward  side  of  a  body-strake  combination,  from  Verle. 
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Cross-Flow  Drag  coefficient 


The  analysis  that  follows  refers  to  and  depends  upon  the  eross-flow  drag  eoefficient  of 
eireular  eylinders  and  two  dimensional  strips,  representing  isolated  strakes.  Thus,  these 
two  basic  shapes  are  being  reviewed  first. 

Circular  Cylinders 

Many  researchers  measured  cross  flow  drag  of  circular  cylinders,  Cdo,  and  several 
authors,  e.  g.  Hoerner,  (Chapters  XV  and  XVI)  Jorgensen,  and  the  British  ESDU 
compiled  available  experimentally  obtained  drag  coefficient  data.  The  latter  will  be  used 
in  the  present  study  because  of  the  extant  of  data  fused  into  it  and  the  accurate 
presentation.  See  Fig  2. 

Plain  Strakes 

The  main  source  for  cross  flow  drag  coefficient  of  a  strake,  Cd|,  is  Lindsey^  who  covered 
Mach  numbers  up  to  0.7.  Hoerner^  compiled  the  available  data.  (Chapter  15-3)  His 
proposed  curve  for  the  subsonic  range  is  presented  in  Fig.  4  that  is  discussed  in  the  next 
section.  The  value  of  Cd|  in  the  subsonic  region  is  about  2.0. 

Circular-Cylinder-Strake  combinations 

As  mentioned  before,  literature  survey  yielded  only  one  study  of  the  cross-flow  drag 
coefficients  of  body-strake  combinations.  Macha  ’  measured  pressure  distributions  on 
plain  circular  cylinders  and  in  combination  with  strakes.  He  integrated  the  pressure 
distributions  to  obtain  drag  coefficients.  His  tests  covered  Mach  number  range  of  0.6  to 
1.05  and  inclusive  strake  height  to  diameter  (=s/d)  ratios  of  1.2,  1.55,  and  2.0.  For  Mach 
numbers  0.6,  0.7,  and  0.8  all  three  stakes  were  tested.  For  Mach  numbers  0.9  and  1.0, 
only  the  narrower  two  were  tested  and  at  Mach  number  of  1.05,  only  one.  Plain  cylinders 
of  various  diameters  were  also  tested  at  Mach  numbers  between  0.6  and  1.2.  The  results 
enabled  him  to  correct  for  blockage  effects.  Since  only  one  diameter  was  tested  with 
strakes,  the  corrections  for  these  cases  were  based  on  the  results  obtained  with  plain 
cylinders. 


Schematic  from  Macha. 


Fig.  2  also  compares  Macha’s  results  for  circular  cylinders  with  the  correlation  of  Ref.  4. 
Good  agreement  is  apparent,  except  for  M=0.9.  This  finding  corroborates  the  quality  of 
these  data. 


Fig.  3,  reproduced  from  Ref  6,  shows  the  dependence  of  cross  flow  drag  coefficients  on 
strake  inclusive  width- to-diame ter  ration,  for  six  Mach  numbers.  It  is  apparent  that 
Macha  included  in  his  charts  plain  cylinder  data  at  s/r=1.0  and  that  his  curves  consider 
these  data.  In  Section  5.5  of  Ref.  6  he  refers  to  the  experimentally  obtained 
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circumferential  pressure  distribution  and  states:  “even  the  smallest  shake  height- to- 
diameter  ratio,  s/d,  of  1.2  the  shake  edge  is  the  dividing  point  between  the  positive  and 
suction  pressures.”  This  eonclusion  is  aceepted  because  the  edges  of  the  shakes 
determine  the  loeation  of  the  separation  points,  which  cause  the  pressure  jump.  It  is 
conjeetured  that  even  very  narrow  shakes  will  force  separation  at  its  edges.  Thus,  the 
present  eurve  fihing  does  not  consider  no-shake  results.  Rather,  straight  lines  were  fitted 
to  the  available  data  and  were  extrapolated  to  s/d=1.0.  Sinee  the  data  are  sparse,  eye-ball 
approximation  was  used.  Since  there  is  only  one  datum  for  Maeh  number  of  1.05,  it  was 
not  ineluded.  The  present  fitted  straight  lines  were  added  to  Fig.  3.  The  available  data  for 
Cd|  at  Mach  numbers  of  0.6  and  0.7  are  also  presented  for  eomparison. 

The  values  of  Cd(|)  at  s/d=1.0  are  named  Cdi  and  those  at  s/d=2.0,  Cd2.  The  present  model 
for  cross  flow  drag  coefficient  of  body-shake  eombinations  is 
Cd^  =Cdi+(Cd2-Cdi)(s/d-1.0), 

=  (2.0-s/d)Cdi+(s/d-1.0)Cd2  (1) 

The  line  for  M=0.9  is  anomalous,  as  is  Cdo  for  that  Maeh  number,  and  was  not  taken  in 
further  considerations.  It  is  also  apparent  that  Cd2  is  close  to  Cd|.  This  implies  that  for 
s/d>2.0  the  cross  flow  drag  coefficient  should  be  about  that  of  a  plain  shake. 

Cdi  and  Cd2  are  presented  in  Fig.  4,  with  Cdo  as  reference.  It  is  apparent  that  they  follow 
the  trends  of  Cdo  the  transonic  region.  The  ratios  Ki=Cdi/Cdo  and  K2=Cd2/Cdo  were 
evaluated  for  the  three  subsonic  Mach  numbers  and  their  averages  are  Ki=1.125  and 
K2=1.34.  The  use  of  these  two  faetors  is  diseussed  in  the  next  seetion. 

Also  shown  in  Fig.  4  is  Cd|  from  Hoerner^.  As  indicated  before,  in  the  subsonie  region, 

Cd2  «Cd|. 
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Jorgensen  ’  estimates  the  cross-flow  drag  coefficient  on  a  winged  (or  straked)  eireular 
body  by  eonsidering  the  normal-foree  acting  on  the  cross  sections,  as  calculated  by  the 
Newtonian  theory. 


t-'d(|)  [CNjjj/CNpjNew  Cdo-  (2) 

With 

[Cn(|,/ CnqInb  =  (3/2)(s/d-l/3),  based  on  d, 

=  (3/2)[1.0-(l/3)(d/s)],  based  on  s.  (3) 

This  relationship  yields  Cdi=Cdo  and  Cd2=1.25Cdo-  Recall  that  the  present  findings  are 
Kl=1.125  and  K2=1.34.  Namely,  the  present  estimate  of  Cd,],  is  larger  than  Jorgensen’s 
by  a  faetor  of  1.125  to  1.07. 

Some  shakes  are  aetually  highly  swept,  low  aspeet  ratio,  A  wings.  Others  feature  sueh 
front  ends,  followed  by  reetangular  rear  parts.  In  the  ease  of  shakes  with  strong  sweep,  it 
is  expected  that  the  added  normal  force,  due  to  separation  along  the  leading  and  the  side 
edges  of  the  shakes,  will  be  determined  by  the  leading-edge  suction  analogy  devised  by 
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Polhamous.^  His  prediction,  at  low  speeds,  is  analogous  to  CopTr.  Leading-edge  suction 
decreases  as  Mach  number  increases  and  diminishes  when  the  leading-edge  becomes 
sonic.  Nevertheless  no  attempt  was  done  to  include  this  effect  in  the  present  analysis. 

Note:  Eqn.  (1)  is  valid  for  s/d<2.0.  For  larger  values  of  s/d,  Cd(|)=Cd2.  However,  large 
values  of  s/d  should  be  considered  wings  rather  than  strakes. 

Adjustment  for  Finite  Bodies 

It  is  known  e.  g.  Welsh  and  Jorgensen,  that  cross  flow  drag  coefficient  of  finite  bodies 
is  smaller  than  that  of  matching  two  dimensional  bodies.  Cross  flow  drag  proportionality 
factor,  defined  as  the  ratio  between  the  cross  flow  drag  coefficients  of  finite  to  infinite 
bodies  was  introduced.  Jorgensen^  processed  his  own  experimentally  obtained  normal- 
force  curves  from  a  previous  study  (Ref.  11)  and  obtained  an  effective  finite  body  cross 
flow  drag  coefficient.  His  results,  from  Fig.  5  of  Ref  3  are  presented  in  Fig.  5  and 
compared  with  Cdo  from  Ref.  4.  It  is  apparent  that  in  the  transonic  region  there  is  a 
noticeable  difference  between  the  cross  flow  drag  coefficients  of  finite  and  infinite 
bodies. 

Based  on  the  available  data,  an  arbitrary  model  was  selected  for  this  work.  It  consists  of 
“sewing”  parts  of  the  available  data: 

Table  1:  Sources  of  the  present  cross  flow  drag  coefficient  for  finite  bodies 


M  range 

source 

0.0 -0.7 

ESDU 

0.7 -0.9 

A  constant  value  connecting  local 
peaks  of  the  two  sources 

0.9-  1.2 

Jorgensen 

1 .2  and  up 

ESDU 

The  present  model  is  presented  in  Fig.  5.  It  features  a  dip  around  M=0.9,  and  a  peak 
around  M=l.l.  No  smoothing  was  introduced  into  the  model. 

Finite  body  models  for  Cdi  and  Cd2  were  obtained  by  multiplying  finite  circular  body 
model  by  Ki  and  K2,  respectfully.  The  subsonic  branch  of  the  latter  was  modified  to 
account  for  the  observed  proximity  of  Cd2  to  Cd|  in  that  range.  The  models  are  depicted  in 
Fig.  6  and  a  tabular  form  is  shown  in  Appendix  A. 
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Cdo  comparison  (2-D) 


Fig,  2:  The  dependence  of  Cdo  on  Mach  number,  from  ESDU  and  Macha, 
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^v\ 


Strcke  height- Cylinder 
Diameter  Ratio-s/d 


Fig,  3:  The  dependence  of  Cd,],  on  stroke  width  to  diameter  ratio,  from  Macha, 
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do> 


3.2 


Cd  of  straked  cylinders  (2-D) 


Fig,  4:  Dependence  of  Cdi  and  Cdi  on  Mach  number. 


Cdo  I  Cdo  comparison 


M 


Fig,  5:  Comparison  between  Jorgensen’s  Cdo  and  Cdo. 
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2.4 


Cqc  (Finite  body) 


- CD1,  Finite  body 

- CD2,  Finite  body 


0.8  I  I  I  I 

0.0  0.4  0.8  1.2  1.6  2.0 

M 


Fig,  6:  Present  cross  flow  drag  coefficient  models. 


14 


Analysis 


Contributions  Buildup 

Following  previous  researchers,  e.  g.  USAF  Datcom^^  and  Hemsch^^,  the  estimation  of 
the  normal-force  and  pitching-moment  coefficients  consider  two  contributions.  The  first 
one  is  that  of  the  inviscid  or  attached  flow,  and  the  second  one  is  due  to  viscous  effects  or 
cross-  flow  separation. 

1 .  The  non  viscous  contributions  are 


Cni  =  V2CNa-sm(2a)  (4. a) 

Cmi=  14  Cm„-sm(2a)  (4.b) 


The  two  stability  derivatives  in  Eqns.  (4)  were  obtained  using  the  Missile  Datcom  code 
(M-Dat,  Ref.  14)  for  zero  angle  of  attack. 


2.  The  present  estimation  of  the  contribution  of  the  separated  cross  flow 
distinguishes  between  body  sections  without  and  with  strakes.  For  the  former, 
namely  for  plain  circular  body,  Cdq  is  taken  into  account,  and  for  the  latter,  Cd(|). 
These  two  drag  coefficients  depend  on  Mc=M  sina.  Thus, 


Where 

f  Cdo  for  plain  body  sections, 

CDe=  ^ 

[  Cd(|)  for  body  sections  with  strakes. 

f  d  for  plain  body  sections, 

B=  ^ 

[  s  for  body  sections  with  strakes. 


(5. a) 
(5.b) 


(6) 

(V) 


Computations 

Test  configurations  were  divided  to  longitudinal  sections  according  to  their  cross 
sectional  shape.  For  triangular  strakes,  analytical  integration  of  Eqns.  (5. a)  and  (5.b)  were 
performed  to  account  for  the  varying  span  and  s/d  ratio.  Eor  details  see  appendix  B.  Eor 
rectangular  strakes,  B=s  and  Eqn.  (1)  becomes  a  constant  along  the  body-strake  sections. 
A  service  code  that  calculates  Me  for  a  given  Mach  number  and  angles  of  attach  and 
provides  values  of  the  three  cross  flow  drag  coefficients  was  used.  The  rest  of  the 
computations  were  done  on  an  excel  worksheets. 
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Validation 


Three  types  of  configurations  were  selected  for  the  validation  of  the  proposed  method: 

1.  Conical  body  and  A  wing  combinations  having  common  apex.  These 
configurations  feature  constant  s/d,  and  common  center  of  pressure  for  both 
components  and  both  contributions,  namely  the  inviscid  and  viscous. 
(Benchmarks  1-2) 

2.  Bodies  and  rectangular,  or  mostly  rectangular,  strakes.  (Benchmarks  3-5) 

3.  Bodies  with  A  strake  combinations.  (Benchmarks  6-7) 


For  the  first  two  benchmarks  Cd,|,  is  constant  along  the  configuration  and  only  depends  on 
Me.  In  these  cases  Eqn.  (5. a)  reduces  to 


CN2  =  (Sp/SR).CDe-sina 


(8) 


The  center  of  pressure  is  at  2/3  the  length. 

Also,  since  these  configurations  are  slender,  generalized  slender  body  theory,  (s-b-t)  e.  g. 
Nielsen,^"^  and  Ashley  and  Landahl,'^  was  used  for  a  second  estimation  of  the  normal- 
force  curve  slope. 

Cn„  =  2(s/d)^[  1  -(d/s)^+(s/d)^]  (9) 


Benchmark  1 

Jorgensen'^,  in  NACA  TN  4045,  tested  a  cone-A-strake  combination  at  Mach  numbers  of 
1.97  and  2.94.  A  schematic  of  his  test  model  is  depicted  in  Fig.  7  and  the  M-Dat 
geometrical  model  is  shown  in  Fig.  8.  This  model  assures  the  nominal  s/d  ratio,  which  is 
needed  for  the  evaluation  of  body-strake  mutual  interference  factors. 

The  fineness  ratio  of  the  conical  body  is  3.67,  and  s/d=1.84. 

The  M-Dat  estimates  of  the  normal-force  curve  slopes  are  given  in  Table  1.  The 
estimated  center  of  pressure  location  is  very  close  to  2/3  body  length,  as  expected. 

Table  2  Normal-force  curve  slope  of  benchmark  1  configuration 


M 

Cno 

M-Dat 

1.97 

4.58 

2.94 

4.33 
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For  s/d=1.84,  Cd^  =  0.16Cdi+0.84Cd2 


Comparisons  between  predictions  and  test  data  are  presented  in  the  two  parts  of  Fig.  9. 
For  M=1.97,  the  agreements  between  the  two  predictions  and  test  data  are  very  good.  For 
M=2.94,  the  two  predictions  are  very  close.  Flowever,  test  data  is  about  linear  (CL  vs.  a) 
while  the  predictions  show  some  nonlinearity,  causing  a  slight  over-estimate  of  the  lift 
coefficient. 


Fig,  7:  Schematic  of  benchmark  configuration  1,  from  Jorgensen, 


Fig,  8:  M-Datcom  geometrical  model  of  benchmark  1, 
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NACA  TN  4045  -  B1W1  -  M=1.97 


0  4  8  12  16  20 

a 

a) 


NACA  TN  4045  -  B1W1  -  M=2.94 


0  4  8  12  16  20 


a 


b) 


Fig,  9:  Comparison  between  predictions  and  test  data  for  benchmark  1, 

a)  M=1.97;  and  b)  M=2.94. 
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Benchmark  2 

Foster'^,  in  NASA  TM-X-167,  tested  a  cone-A-strake  combination  at  Mach  numbers  of 
1.41  and  2.01.  A  schematic  of  his  test  model  is  depicted  in  Fig.  10.  The  conical  body  has 
fineness  ratio  of  4.0,  and  s/d=2.0,  yielding  Sp/So=10.19.  The  M-Dat  geometrical  model  is 
shown  in  Fig.  11.  The  estimates  of  the  normal-force  curve  slopes  are  given  in  Table  2. 

Table  3  Normal-force  curve  slope  of  benchmark  2configuration 


M 

CNa 

M-dat 

s-b-t 

1.41 

6.07 

6.5 

2.01 

5.95 

6.5 

Since  s/d=2.0,  Cd,|,  =  Cd2- 

Comparisons  between  analysis  and  test  data  are  shown  in  Fig.  12.  For  M=1.41,  M-Dat 
predicts  Cn  which  is  up  to  10.0%  lower  than  the  data.  The  present  estimates  are  up  to 
15.0%  lower.  For  M=2.01  the  two  predictions  are  in  very  good  agreement  with  the  data. 


Fig,  10:  Schematic  of  benchmark  2  configuration,  from  Foster, 


NASA  TM-X-167 


Fig,  11:  M-Datcom  geometrical  model  of  benchmark  2, 
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NASATM-X-167  Models  M=1.41 


a) 


b) 

Fig,  12:  Comparison  between  predictions  and  test  data  for  benchmark  2, 

a)  M=1.41;  and  b)  M=2.01. 
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Benchmark  3 

Reference  1 1  by  Jorgensen  contains  the  database  that  was  used  by  him  to  devise  the 
effect  of  finite  body  length  on  the  cross  flow  drag  coefficient.  In  Ref.  19  he  tested  the 
same  bodies  in  combination  with  narrow  strakes.  Configuration  Nl-Cl-S  was  selected 
and  is  presented  in  Fig.  13.  The  geometry  of  the  input  model  is  shown  in  Fig.  14.  The 
chord  of  the  strake  matches  the  length  of  the  cylindrical  afterbody.  Span  to  diameter  ratio 
is  1.2.  Reference  point  for  center  of  pressure  location  is  body  base. 

The  stability  derivatives  of  the  body  alone  and  of  the  configuration,  as  obtained  by  the  M- 
Dat  code  are  given  in  table  4. 

Table  4  Stability  derivatives  of  benchmark  3  configuration 


Body  alone 

Body-strake 

M 

Cnh 

Cma 

Cnq 

Cma 

0.9 

2.49 

18.07 

2.94 

20.40 

1.2 

2.61 

19.57 

2.97 

21.10 

1.5 

2.48 

20.97 

6.80 

39.22 

2.0 

2.75 

23.27 

5.52 

24.98 

It  is  apparent  that  the  normal-force  curve  slope  of  the  configuration  is  much  larger  at  the 
two  highest  Mach  numbers  than  those  of  the  low  ones.  This  is  attributed  to  a  jump  in  the 
contribution  of  the  fins  alone  between  Mach  number  1.2  and  1.5.  For  the  transonic  Mach 
numbers,  the  normal-force  curve  slope  of  the  isolated  fins  is  close  to  those  predicted  by 
slender  wing  theory.  In  the  supersonic  range,  the  values  obtained  by  the  code  are  an  order 
of  magnitude  larger  than  those  expected  by  the  theory.  (The  aspect  ratio  of  a  pair  of 
strakes  is  about  0.03.) 

Comparisons  between  predictions  and  test  data  are  presented  in  the  first  four  parts  of  Fig. 
15,  for  Mach  numbers  of  0.9,  1.2,  1.5  and  2.0.  At  the  two  low  Mach  numbers,  the  present 
method  predicts  well  the  normal-force  coefficient,  including  the  local  trends,  to  angles  of 
attack  of  50  deg.  The  estimate  of  the  M-Dat  is  about  15%  lower  than  the  data.  At  the 
higher  Mach  numbers,  M-Dat  and  present  predictions  of  the  normal-force  coefficient  at 
small  and  moderate  angles  of  attack  are  higher  than  test  data.  At  large  angles  of  attack, 
the  present  prediction  departs  from  the  test  data  and  remains  higher  than  the  data.  The 
departure  starts  at  angles  of  attack  where  the  cross  flow  Mach  number  is  about  1.0.  At 
this  Mach  number  the  cross  flow  drag  coefficients  feature  a  second  hump.  Generally,  the 
center  of  pressure  location  is  well  predicted  by  both  methods,  except  for  the  dip  which  is 
characteristic  for  the  two  low  Mach  numbers. 

The  difference  in  the  initial  slopes  of  the  normal-force  curves,  which  was  found  at  M=1.5 
and  M=2.0,  develops  noticeable  differences  at  high  angles  of  attack  as  well.  An 
independent  estimate  of  the  slope  was  sought.  Literature  search  for  reliable  data  for  very 
low  aspect  ratio  rectangular  wings  did  not  yield  result.  Thus,  a  hybrid  estimate  was  used. 
The  contribution  of  the  body  alone  was  obtained  from  the  M-Dat  code.  The  contribution 
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of  the  strake  unit  (strake  and  mutual  interference  with  the  body  was  estimated  using 
generalized  s-b-t. 


Cn„(B-S) 

=  2-(b/d)^-[l-(d/b)^+(d/b)^] 

(lO.a) 

Cn„  (B) 

=  2.0 

(lO.b) 

Cn„(SU) 

=  CNa(B-S)  -  CNa  (B) 

Cn„(SU) 

=  2-(b/d)^-[l-2(d/b)^+(d/b)^] 

Cn„(SU) 

=  2-(b/d)^-[l-(d/b)^]^ 

(11) 

For  b/d  =  1 .2,  Cnk  (SU)  =  0.269.  Thus,  for  the  present  case 


Cn„(B-S)  =  Cn„(B)  + 0.269. 


(12) 


The  hybrid  estimate  of  Cnh  is  presented  in  Table  5.  For  the  lower  Mach  numbers,  the 
agreement  between  the  two  estimates  is  good.  For  the  higher  Mach  numbers,  the  hybrid 
estimates  are  smaller  than  predicted  by  the  M-Dat  code.  Parts  (e)  and  (f)  of  Fig.  15  show 
repeated  comparisons  between  test  data  and  the  present  method  with  Cnk  from  both  M- 
Dat,  (Present- 1)  and  revised  Cnh.  (Present-2)  It  is  apparent  that  the  estimates  that  are 
based  on  the  hybrid  approach  improves  the  agreement  with  test  data  at  small  and 
moderate  angles  of  attack. 

Table  5  Comparison  of  normal-force  curve  slopes 


M 

B,  M-Dat 

S,  M-Dat 

B-S,  M-Dat 

B-S,  Hybr, 

0.9 

2.49 

0.19 

2.94 

2.76 

1.2 

2.61 

0.09 

2.97 

2.88 

1.5 

2.48 

1.27 

6.80 

2.75 

2.0 

2.75 

0.82 

5.52 

3.02 

At  larger  angles  of  attack,  the  revised  Cn  curves  still  deviate  from  the  test  data.  As  noted 
before,  the  deviation  occurs  around  Mc=1.0.  This  may  imply  that  the  cross  flow  drag 
coefficient  of  the  present  model  is  a  little  too  high  at  low  supersonic  Mach  numbers. 
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Fig,  13:  Schematic  of  benchmark  3  configuration,  from  Jorgensen, 


TMX-3130  N1-C1-S 

TOP  VIEW  PROFILE  OF  F1  =  HEXAGONA 


Fig,  14:  M-Datcom  geometrical  model  for  benchmark  3, 
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20 


NASATM-X-3130,  M=0.9 


NASATM-X-3130,  M=0.9 


Fig,  15:  Comparison  between  predictions  and  test  data  for  benchmark  3, 

a)  M=0.9. 
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0  10  20  30  40  50  60 
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NASATM-X-3130,  M=1.2 


Fig,  15:  Comparison  between  predictions  and  test  data  for  benchmark  3, 

b)  M=1.2. 
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0  10  20  30  40  50  60 
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NASATM-X-3130,  M=1.5 


Fig,  15:  Comparison  between  predictions  and  test  data  for  benchmark  3, 

c)  M=1.5. 
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0  10  20  30  40  50  60 
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NASATM-X-3130,  M=2.0 


Fig,  15:  Comparison  between  predictions  and  test  data  for  benchmark  3, 

d)  M=2.0. 
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0  10  20  30  40  50  60 

a 

Fig,  15:  Repeated  comparison  between  predictions  and  test  data,  benchmark  3, 

e)  M=1.5. 


0  10  20  30  40  50  60 


a 


Fig,  15:  Repeated  comparison  between  predictions  and  test  data,  benchmark  3, 

f)  M=2,0, 
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Benchmark  4 

From  the  many  configurations  that  were  tested  by  Robinson,  in  NACA  RM  L57D19, 
the  three  that  feature  rectangular  shakes  were  selected.  Mach  number  of  his  tests  was 
2.01.  Schematics  of  the  body  and  the  shakes  are  depicted  in  the  two  parts  of  Fig.  17.  The 
body  consists  of  a  tangent  ogive  nose  with  a  fineness  ratio  of  3.5  and  a  cylindrical 
afterbody  with  fineness  ratio  of  6.5.  All  shakes  have  chord  to  diameter  ratio  of  4.333. 
Their  span-to-diameter  ratios  are  given  in  Table  4. 

Table  5  Strake  span  to  diameter  ratios  for  benchmark  4  configurations 


Configuration 

b/d 

F7-W6-A1 

4/3=1.333 

F7-W5-A1 

5/3=1.667 

F7-W4-A1 

7/3=2.333 

A  typical  M-Dat  geometrical  model  is  shown  in  Fig.  18. 

Comparison  between  test  data  and  predictions  for  the  body  alone  are  shown  in  Fig.  19.  At 
small  angles  of  attack  the  two  methods  overestimate  the  test  data.  At  angles  of  aback 
larger  than  14.0,  the  present  prediction  of  the  normal-force  coefficient  is  closer  to  the  data 
than  that  of  the  M-Dat.  The  two  estimates  of  the  pitching-moment  curve  are  in  very  good 
agreement,  however,  their  nonlinearity  are  not  as  pronounced  as  that  of  the  test  data.  This 
accounts  for  a  difference  of  0.35  diameters  in  the  center  of  pressure  location. 

Preliminary  comparison  between  predictions  and  test  data  revealed  that  the  predicted 
normal-force  curve  slope  of  configuration  F7-W4-A1  is  about  28%  lower  than 
experimentally  observed.  As  a  result,  the  agreement  between  prediction  and  test  data  is 
poor  at  all  angles  of  attack.  Thus,  this  configuration  is  not  included  in  the  discussion  that 
follows. 

Fig.  20  contains  comparisons  between  estimates  and  test  data  for  the  configurations  with 
the  two  narrow  shakes.  The  two  predictions  are  very  close.  The  normal-force  coefficient 
for  configuration  with  shakes  W5  and  W6  are  very  close  to  test  data.  The  pitching- 
moment  curves  show  same  trend  as  body  alone,  with  about  the  same  difference  in  the 
center  of  pressure  location. 
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F^A, 


Rectongulof  wing  series 

Fig,  16:  Schematic  of  the  modular  hody-strake  model,  from  Rohinson. 


RM  A57D19  F7-W5-A1 


Fig,  17:  M-Datcom  input  model  for  configuration  F7-W5-A1, 


0  0.2  0.4  0.6  0.8  1 


Me 


Fig,  18:  Cross  flow  drag  coefficient  for  the  configurations  of  benchmark  4, 
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F7-W0-A1  -  Body  alone 


a) 


F7-W0-A1  -  Body  alone 


b) 

Fig,  19:  Comparison  between  predictions  and  test  data  for  the  body  of  benchmark  4, 
a)  normal-force  coefficient,  and  b)  pitching-moment  coefficient. 
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NACARM  L57D19,  F7-W*-A1,  M=2.01 


Fig,  20:  Comparison  between  predictions  and  test  data  for  the  body-strake 
combinations  of  benchmark  4,  a)  normal-force  coefficient. 
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NACA  RM  L57D19,  F7-W*-A1,  M=2.01 


Fig,  20:  Comparison  between  predictions  and  test  data  for  the  body-strake 
combinations  of  benchmark  4,  b)  pitching-moment  coefficient. 
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Benchmark  5 

Simpson  and  Birch  tested  a  fineness  ratio  19.0  body  with  three  very  low  aspeet  ratio 
wings.  Wing  22,  which  is  actually  a  strake  with  b/d=1.6  was  selected  for  this  work.  A 
schematic  of  the  configuration  and  the  computational  model  are  depleted  in  Figs.  21  and 
22,  respectively. 


BllA 


- 

1 

3D 

16D 

B11AW22A3 

<1 

-  1 

3D 

Fig,  21:  Schematic  of  benchmark  5  configuration,  from  Simpson  and  Birch, 


Fig,  22:  Computational  model  for  benchmark  5  configuration. 

Comparisons  between  calculations  and  test  data  are  presented  in  Fig.  23.  For  the  body 
alone,  the  two  predietions  are  very  elose  and  in  good  agreement  with  the  test  data.  In  the 
case  of  the  combination,  the  present  estimate  of  the  normal-force  coefficient  is  in  very 
good  agreement  with  the  data,  while  the  M-Dat  over-prediets  it.  The  predicted  center  of 
pressure  loeation  is  more  forward  than  experimentally  observed,  with  small  differences 
between  the  two  methods.  The  maximum  gap  is  about  1.5  body  diameters  at  moderate 
angles  of  attack  and  slightly  better  at  higher  angles  of  attack. 
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AIAA-2001-2410 


0  4  8  12  16  20  24 

a 


a) 


AIAA-2001-2410 


b) 

Fig,  23:  Comparison  between  predictions  and  test  data  for  benchmark  5, 
a)  Normal-force  coefficient,  and  b)  center  of  pressure  location. 
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Benchmark  6 

Robinson  investigated  at  a  Mach  number  of  2.01  a  fineness  ratio  10.0  body  with  several 
wings  and  controls.  Among  others,  he  tested  a  5.0  deg  delta  strake  having  b/d=  2.07.  A 
schematic  of  that  configuration  is  shown  in  Fig.  24. 

Since  b/d>2.0,  the  strake  was  divided,  lengthwise,  into  two  part.  The  front  one  ending  at 
body  station  where  b/d=2.0  was  analyzed  according  to  the  expression  developed  in 
Appendix  A.  For  the  rear  part  Cd2  was  considered  as  cross  flow  drag  coefficient. 

Fig.  26  is  a  comparison  between  predictions  and  test  data  for  this  configuration.  The  two 
predictions  of  the  normal-force  coefficient  are  close  and  match  the  data,  except  for  the 
present  method  that  slightly  overestimates  the  data  past  angle  of  attack  of  22.0  deg.  There 
is  a  slight  difference  between  the  two  estimates  of  the  pitching-moment  curve  and  the 
data.  The  present  method  predicts  stronger  nonlinearity  in  this  curve  and  better  predicts 
the  data  at  high  angles  of  attack. 


Fig,  24:  Schematic  of  body  with  5“  strake,  from  Robinson. 


Fig,  25:  M-Datcom  geometrical  model  for  Benchmark  6, 
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NACA  RM  L58A21  5°  delta  strake 


a) 


NACA  RM  L58A21  5°  delta  strake 


a 

b) 

Fig,  26:  Comparison  between  predictions  and  test  data  for  benchmark  6, 
a)  normal-force  coefficient,  and  b)pitching-moment  coefficient. 


37 


Benchmark  7 

Spearman  and  Robinson  tested  a  missile  composed  of  a  fineness  ratio  10.0  body  with 
low  aspect  ratio  wing  and  aft  control.  As  can  be  seen  form  Fig.  27,  the  control  surfaces 
are  very  close  to  the  trailing-edge  of  the  wing,  thus  were  treated  as  one  unit.  The 
geometrical  model  is  depicted  in  Fig.  28.  This  benchmark  covers  Mach  numbers  of  2.01, 
4.65  and  6.8. 

Comparisons  between  predictions  and  test  data  are  presented  in  the  three  parts  of  Fig.  29. 
For  M=2.01,  the  two  predictions  are  close  and  in  good  agreement  with  test  data.  The 
difference  in  Cm  is  equivalent  to  a  fraction  of  body  diameter  in  the  center  of  pressure 
location.  For  the  two  high  Mach  numbers,  the  analysis  overestimates  Cn  and  the  absolute 
value  of  Cm.  In  this  case  too,  the  present  analysis  deviates  from  the  trend  of  the  data  as 
Mc=1.0  and  higher. 


Fig,  27:  Schematic  of  benchmark  7  configuration,  from  Spearman  and  Robinson, 


TM  X-46 


Fig,  28:  M-Datcom  input  model  for  benchmark  7, 
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10 


NASATM-X-46,  M=2.01 


NASATM-X-46,  M=2.01 


Fig,  29:  Comparison  between  predictions  and  test  data  for  benchmark  7, 

a)  M=2.01. 
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NASATM-X-46,  M=4.65 


Fig,  29:  Comparison  between  predictions  and  test  data  for  benchmark  7, 

b)  M=4.65. 
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Fig,  29:  Comparison  between  predictions  and  test  data  for  benchmark  7, 

c)  M=6.8. 
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Summary  and  Remarks 


1.  The  longitudinal  characteristics  of  body-strake  configurations  were  evaluated 
using  component  buildup  approach  that  features  a  unified  approach  for  the 
estimation  of  the  vortex  induced  contributions  of  both  components.  The  new 
approach  uses  one  cross  flow  drag  coefficient  for  the  part  of  the  body  on  which 
the  strakes  are  mounted.  The  basis  for  the  drag  model,  are  empirical  data  of 
Macha  who  tested  straked  cylinders  in  two-dimensional  transonic  flow. 
Corrections  were  made  to  account  for  finite  body  effect  and  for  the  cross  flow 
drag  coefficient  of  a  plain  shake  at  subsonic  Mach  numbers. 

2.  A  component  buildup  method  was  formulated,  that  combines  non  viscous  and 
viscous  contributions  to  the  normal-force  and  pitching-moment  coefficients.  For 
the  first  contribution,  the  stability  derivatives  were  obtained  using  the  Missile 
Datcom  code.  A  unified  analysis  was  used  for  the  other  contribution,  using  the 
present  cross  flow  drag  coefficient  model. 

3.  A  pre-processing  code  was  written  to  obtain  the  cross  flow  drag  coefficient  of 
each  case  as  a  function  of  cross  flow  Mach  number  and  span  to  diameter  ratio. 
Analytical  integration  was  done  for  delta  strakes  where  span  to  diameter  ratio 
varies  along  the  chord. 

4.  The  new  method  was  evaluated  by  comparing  predictions  with  test  data 
assembled  from  the  literature.  Predictions  by  the  Missile  Datcom  code  were  also 
included  for  additional  comparison.  In  most  cases,  the  agreements  between  the 
Missile  Datcom  and  the  present  method  are  good  or  very  good.  In  the  case  of  two 
configurations  that  feature  rectangular  strakes,  the  present  analysis  provides 
estimates  that  are  much  closer  to  test  data  than  those  of  the  Missile  Datcom  code. 

5.  In  one  case,  that  features  a  very  low  aspect  ratio  shake,  the  Missile  Datcom  code 
provided  uncertain  estimates  for  the  normal-force  curve  slope  of  the  combination, 
which  is  attributed  to  that  of  shake  alone.  A  hybrid  approach  that  uses  generalized 
slender  body  theory  and  Missile  Datcom  characteristics  for  body  alone, 
considerably  improved  agreement  between  estimates  and  test  data. 

6.  The  present  method  gives  good  estimates  as  long  as  the  cross  flow  Mach  number 
is  smaller  than  0.95.  At  higher  values,  the  predictions  deviate  from  the  data  and 
the  predicted  normal-force  coefficient  is  higher  than  experimentally  obtained. 
This  may  imply  that  the  present  cross  flow  drag  coefficient  is  a  little  too  high  at 
low  supersonic  Mach  numbers.  It  is  suggested  to  use  additional  test  data  to  adjust 
the  present  drag  model. 

7.  There  is  a  need  for  additional  test  data  for  body-strake  combinations  at  tri-sonic 
Mach  numbers. 


42 


List  of  References 


1.  Werle,  H.,  “Tourbillons  de  Corps  Fusseles  aux  Incidence  Eleve,”  I’Aeronautique 
et  I’Astronautique,  JNe  79,  June  1979. 

2.  Hoemer,  S.  F.,  “Fluid  Dynamics  Drag,”  Author’s  Edition,  Midland  Park,  NJ, 
1965. 

3.  Jorgensen,  F.  H.,  “Prediction  of  Static  Aerodynamic  Characteristics  for  Slender 
Bodies  Alone  and  With  Fifting  Surfaces  to  Very  High  Angles  of  Attack,”  NASA 
TRR-474,  Sep.  1977. 

4.  Anon.,  “Drag  of  Circular  Cylinders  Normal  to  a  Flat  Plate  With  a  Turbulent 
Boundary  Fayer  for  Mach  Numbers  Up  to  3,”  ESDU  International,  Item  83025, 
Fondon,  UK,  1983. 

5.  Findsey,  W.  F.,  “Drag  of  Cylinders  of  Simple  Shapes,”  NACA  Rep.  619,  Oct. 
1937. 

6.  Macha,  J.  M.,  “A  Wind  Tunnel  Investigation  of  Circular  and  Staked  Cylinders  in 
Transonic  Cross  Flow,”  TAMU  Rep.  3318-76-01,  Texas  A&M  Research 
Foundation,  College  Station,  TX,  Nov.  1976. 

7.  Macha,  J.  M.,  “Drag  of  Circular  Cylinders  at  Transonic  Mach  Numbers,”  JA,  Vol. 
14,  No.  6,  June  1977,  pp.605-607. 

8.  Jorgensen,  F.  H.,  “A  Method  for  Estimating  Static  Aerodynamic  Characteristics 
for  Slender  Bodies  circular  and  Noncircular  Cross  Section  Alone  and  With  Fifting 
Surfaces  at  Angles  of  Attack  from  0°  to  90“,”  NASA  TN  D-7228,  Apr.  1973. 

9.  Polhamus,  E.  C.,  “A  Concept  of  the  Vortex  Fift  of  Sharp-Edge  Delta  Wings 
Based  on  a  Feading-Edge-Suction  Analogy,”  NASA  TN  D-3767,  Dec.  1966. 

10.  Welsh,  C.  J.,  “The  Drag  of  Finite-Fength  Cylinders  Determined  from  Flight  Tests 
at  High  Reynolds  Numbers  for  a  Mach  Number  Range  from  0.4  to  1.3,”  NACA 
TN2941,  1952. 

11.  Jorgensen,  F.  H.,  and  Nelson,  E.  R.,  “Experimental  Aerodynamic  Characteristics 
for  a  Cylindrical  Body  of  Revolution  with  Various  Noses  at  Angles  of  Attack 
from  0°  to  58“  and  Mach  Numbers  from  0.6  to  2.0,”  NASA  TM  X-3128,  Dec. 
1974. 


43 


12.  Hoak,  D.  E.,  (Project  Engineer)  and  Ellison,  D.  E.,  (Chief  Scientist)  “USAE 
Stability  and  Control  Dactom,”  Elight  Control  Division,  AE  Elight  Dynamics 
laboratory,  (now  AERL)  Wright-Patterson  AE  Base,  OH.  1968. 

13.  Hemsch,  M.  J.,  “Component  Build-Up  Method  for  Engineering  Analysis  of 
Missiles  at  Eow-to-High  Angles  of  Attack.”  Chapter  4  in  “Tactical  Missiles 
Aerodynamics:  Prediction  Methods,”  Progress  in  Astronautics  and  Aeronautics, 
Vol.  142,  AlAA  Publication,  1992. 

14.  Blake,  W.  B.,  “Missile  Datcom  User’s  Manual  -  1997  Eortran  90  Version,” 
USAE  Research  Eaboratory,  Report  AERL-VA-WP-TR-1988-3009,  Wright- 
Patterson  AEB,  OH,  Eeb.  1998. 

15.  Nielsen,  J.,  “Missile  Aerodynamics,”  NEAR  Inc.,  Mountain  View,  CA,  1988. 

16.  Ashley,  H.,  and  Eandahl,  M.,  “Aerodynamics  of  Wings  and  Bodies,”  Addison 
Wesley  Eongman,  Reading,  MA,  1965. 

17.  Jorgensen,  E.  H.,  “Elliptic  Cones  Alone  and  with  Wings  at  Supersonic  Speeds,” 
NACATN4045,  1957. 

18.  Eoster,  G.  V.,  “Static  Stability  Characteristics  of  a  Series  of  Hypersonic  Boost- 
Glide  Configurations  at  Mach  Numbers  1.41  and  2.01,  “NASA  TM  X-167,  Aug. 
1957. 

19.  Jorgensen,  E.  H.,  “Experimental  Aerodynamic  Characteristics  for  a  Cylindrical 
Body  of  Revolution  With  Side  Strakes  and  Various  Noses  at  Angles  of  Attack 
from  0°  to  58“  and  Mach  Numbers  from  0.6  to  2.0,”  NASA  TM  X-3130,  March 
1975. 

20.  Robinson,  R.  B.,  “Aerodynamic  Characteristics  of  Missile  Configurations  with 
Wings  of  Eow  Aspect  Ratio  for  various  Combinations  of  Eorebodies,  Afterbodies, 
and  Nose  Shapes  for  Combined  Angles  of  Attack  and  Sideslip  at  a  Mach  Number 
of2.01,”NACARM  A57D19,  Jun.  1957. 

21.  Simpson,  G.  M.,  “Aerodynamic  Characteristics  of  Missiles  having  very  Eow 
Aspect  Ratio  Wings,”  AIAA2001-2410,  June  2001. 

22.  Robinson,  R.  B.,  “Wind-Tunnel  Investigation  at  a  Mach  Number  of  2.01  of  the 
Aerodynamic  Characteristics  in  Combined  Angle  of  Attack  and  Sideslip  of 


44 


Several  Hypersonic  Missile  Configurations  With  Various  Canard  Controls,” 
NACA  RM  L58A21,  March  1958. 

23.  Spearman  M.  L.,  and  Robinson,  R.  B.,  “Longitudinal  Stability  and  Control 
Characteristics  at  Mach  Numbers  of  2.01,  4.65,  and  6.8  of  Two  Hypersonic 
Missile  Configurations,”  NASA  TM  X-46,  Sep.  1959. 


45 


Appendix  A 


Tabulated  cross  flow  drag  coefficients  for  finite  bodies 

Present  model  for  finite  length  bodies  and  parameters  for  the  model  for  body-strake 
combinations.  (See  Eqn.  1) 

Cd2  was  modified  in  the  subsonic  range  to  account  for  shake  alone  cross  flow  data. 


Me 

Cdo 

Cdi 

Cd2 

0.00 

1.202 

1.3523 

1.9900 

0.04 

1.202 

1.3523 

1.9900 

0.08 

1.202 

1.3523 

1.9900 

0.12 

1.204 

1.3545 

1.9900 

0.16 

1.207 

1.3579 

1.9900 

0.20 

1.211 

1.3624 

1.9900 

0.24 

1.219 

1.3714 

1.9900 

0.28 

1.226 

1.3793 

1.9900 

0.32 

1.241 

1.3961 

1.9900 

0.36 

1.260 

1.4175 

1.9900 

0.40 

1.280 

1.4400 

1.9900 

0.44 

1.304 

1.4670 

1.9980 

0.48 

1.338 

1.5053 

2.0060 

0.52 

1.371 

1.5424 

2.0160 

0.56 

1.416 

1.5930 

2.0280 

0.60 

1.470 

1.6538 

2.0400 

0.64 

1.528 

1.7190 

2.0520 

0.68 

1.576 

1.7730 

2.1118 

0.70 

1.580 

1.7775 

2.1172 

0.72 

1.580 

1.7775 

2.1172 

0.74 

1.580 

1.7775 

2.1172 

0.76 

1.580 

1.7775 

2.1172 

0.78 

1.580 

1.7775 

2.1172 

0.80 

1.580 

1.7775 

2.1172 

0.82 

1.580 

1.7775 

2.1172 

0.84 

1.580 

1.7775 

2.1172 

0.86 

1.580 

1.7775 

2.1172 

0.88 

1.580 

1.7775 

2.1172 

0.90 

1.530 

1.7213 

2.0502 

0.92 

1.523 

1.7134 

2.0408 

0.93 

1.521 

1.7111 

2.0381 

0.94 

1.520 

1.7100 

2.0368 

0.95 

1.520 

1.7100 

2.0368 

0.96 

1.520 

1.7100 

2.0368 

0.98 

1.525 

1.7156 

2.0435 
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Me 

Cdo 

Cdi 

Cd2 

1.00 

1.530 

1.7213 

2.0502 

1.02 

1.545 

1.7381 

2.0703 

1.04 

1.565 

1.7606 

2.0971 

1.06 

1.590 

1.7888 

2.1306 

1.08 

1.635 

1.8394 

2.1909 

1.10 

1.645 

1.8506 

2.2043 

1.12 

1.650 

1.8563 

2.2110 

1.14 

1.650 

1.8563 

2.2110 

1.16 

1.640 

1.8450 

2.1976 

1.18 

1.630 

1.8338 

2.1842 

1.20 

1.610 

1.8113 

2.1574 

1.24 

1.576 

1.7730 

2.1118 

1.28 

1.542 

1.7348 

2.0663 

1.32 

1.516 

1.7055 

2.0314 

1.36 

1.492 

1.6785 

1.9993 

1.40 

1.470 

1.6538 

1.9698 

1.44 

1.450 

1.6313 

1.9430 

1.48 

1.436 

1.6155 

1.9242 

1.52 

1.419 

1.5964 

1.9015 

1.56 

1.408 

1.5840 

1.8867 

1.60 

1.400 

1.5750 

1.8760 

1.70 

1.376 

1.5480 

1.8438 

1.80 

1.359 

1.5289 

1.8211 

1.90 

1.343 

1.5109 

1.7996 

2.00 

1.327 

1.4929 

1.7782 
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Appendix  B 


The  Non  Linear  Normal-Force  Coefficient  of  Body-A-Strake  Combination 


Objective 

The  nonlinear  eontributions  to  the  normal-foree  and  pitehing-moment  eoeffieients  as 
formulated  in  Chapter  Analysis  are 


Cn2  -  o  J  Cdc  B  dx  •  sin^a  /Sr 

'f  2 

Cm2  -  o  J  Cdc  B  (x-xr)  dx  '  sln  a  /(dRSR) 

Where 

f  Cdo  for  plain  body  sections, 

Cdc=  ^ 

[  Cd(|)  for  body  sections  with  strakes. 

f  d  for  plain  body  sections, 

B=  ^ 

[  s  for  body  sections  with  strakes. 

And 

Cd(|)  =  (2  -  s/d)-CDi  +  (s/d  -  1)'Cd2 


(B.l.a) 

(B.l.b) 


(B.2) 


(B.3) 

(B.3) 


The  integration  of  Eqns.  (B.l.a)  and  (B.l.b)  will  be  only  performed  along  the  chord.  For 
simplicity  the  range  of  integration  will  be  from  o  to  c  and  the  reference  point  will  be  the 
apex  of  the  strake.  The  actual  location  of  the  strake  will  be  taken  into  account  during 
application. 

Define 


In  -  o  J  Cdc  B  dx 

Im  ~  o  j  Cdc  B  X  dx 


(B.4.a) 

(B.4.b) 


The  local  span  is  related  to  the  maximum  span  by 
s  =  d  +  (b-d)-x/c 


(B.5) 
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Define 


a  = 

(b-d)/c 

(B.6.a) 

A  = 

(b-d)/d 

(B.6.b) 

Hence 

s  = 

d  +  a-x 

(B.7.a) 

s/d  = 

1  +  a-x/d 

(B.7.b) 

Introducing  Eqn.  (B.7.b)  into  Eqn.  (B.3)  and  Cdc  B  =  Cd(|)  s,  yields 

Cd(|) 

(1-  a-x/d) -Cdi  +  (a'x/d)  •Cd2 

(B.8) 

Eqns.  (B.4)  become 

In  - 

0  J  [(1 .0-a-x/d)  •  Cdi  +  (a-x/d)  •  Cd2]  (d  +  a-x)  dx 

C  - 

(B.9.a) 

Im  = 

0  J  [(1 .0-a-x/d)  •  Cdi  +  (a-x/d)  •  Cd2]  (d  +  a-x)  x  dx 

(B.9.b) 

Expanding  and  integrating  these  expressions  and  using  definition  (B.6.b)  gives 

In  - 

(c-d)  -  [( 1  -  CbA")  •  Cdi  +  a  •  (1/2  +  /2A)  •  Cd2  ] 

(B.lO.a) 

Im  = 

(c"-d)  -  [(/2  -  %A^)  •  Cdi  +  A  •  (1/3  +  /4A)  •  Cd2  ] 

(B.lO.b) 

Since  Sr  = 

(7t:/4)  •  d  ,  this  gives 

Cn2  = 

{AIk)  -  (c/d)-lN-sinV 

(B.ll.a) 

Cm2  ~ 

{Ain)  -  (c/d)^-Im-sin^a 

(B.ll.b) 
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Appendix  C: 

Planform  Area  and  Center  of  Area  for  Tangent  Ogives 

Geometry: 

Define: 

f=  fn/d  (nose  fineness  ratio) 

Fi= 

F2  =  f  ^  -  % 


From  the  schematic,  it  is  apparent  that 

sin(0o)  =  f  n  /  R 
cos(0o)  =  (R-d  /  2)  /  R 

2  2 

using  sin  0  +  cos  0=1,  one  obtains 
R=  Fi-d 


The  Ogive  Projection  Area 

Ss  =  V2  •  00  •  R^  (Circular  Sector) 

Sa=  V2  •  In  ■  (R-d/2)  (Triangle  defined  by  ogive  axis  and  confining  radii) 

Giving 

Sp=  00 -R^-fn- (R-d/2) 

In  the  parametric  form 

Sp=  d"-(Fi"-0o-f-F2) 

The  Ogive  Moment  of  Area  relative  to  the  base 

The  aera  moment  relative  to  the  base  of  the  ogive  is 

Mo  =  o  I  y(x)  •  X  dx 

X  =  R  sin0 

y  =  R  COS0  -  (R-d/2) 

dx  =  R  COS0  d0 

00  j* 

Mo  =0  J  {R  COS0  -  (R-d/2)}  •  R  cos0  d0 
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Integration  yields 

Mo  =  '/2  •  (R-d/2)  •  (eos(20o)  -  1)  -  %  R^  •  (cos^Bo  -  1) 

Substituting  2eos  0o-l  for  oos20o,  and  replaeing  R  with  the  parametrie  form, 

Mo  =  d^  •  {  F2  •  (F2"  -  Fi')  -  %  •  (F2^  -  Fi^)} 

Normalized  form: 

The  projeetion  area  is  normalized  by  the  cireumscribing  reetangle.  The  eenter  of  area  is 
normalized  by  the  length  of  the  ogive,  and  is  measured  with  respect  to  the  base  of  the 
ogive. 

Sp/(f„-d)  =0o-FiVf-F2 

Xc  /  fn  =  Mo/(Sp  •  f„)  =  {F2  •  (F2'  -  Fi')  -  %  •  (F2'  -  Fi')}  /  {f  •  (Fi^  •  00  -  f  •  F2)} 
The  normalized  forms  are  presented  in  Fig.  Cl. 


Tangent  Ogive  Geometry 


f 


Fig,  Cl:  The  normalized  projection  area  and  center  of  area  of  a  tangent  ogive. 

The  asymptotic  values  are  2/3  for  the  area  ratio  and  3/8  for  the  normalized  area  center. 


51 


